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1^ (57) Abstract: Estimating a frequency of a sampled cardiac rhythm signal and classifying the rhythm. The received signal is sampled 
^ and transformed into a curvature series. A lobe in the curvature series corresponds to a characteristic point in the sampled series. 
O Characteristic points arc selected based on a time of a lobe in the curvature series and, in one embodiment, an amplitude of the signal 
^ at the time of the lobe. A frequency of the sampled series is estimated by autocorrelating a function of the series of the characteristic 
Q points. In one embodiment, the function is a time difference function. The rhythm is classified by plotting the timewise proximity 
of characteristic points derived from an atrial signal with characteristic points derived from a ventricular signal. Regions of the plot 
are associated with a particular rhythm and the grouping of the data corresponds to the classification. 
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TACHYARRHYTHMIA DETECTION AND DISCRIMINATION 
BASED ON CURVATURE PARAMETERS 

5 

Related Applications 

This document is related to co-pending and commonly assigned U.S. 

Patent Application Serial No. 09/703,269, entitled "CURVATURE BASED 

METHOD FOR SELECTING FEATURES FROM AN 
1 0 ELECTROPHYSIOLOGICAL SIGNALS FOR PURPOSE OF COMPLEX 

IDENTIFICATION AND CLASSIFICATION," inventor Sweeney et al, filed 

on October 31, 2000, the specification of which is hereby incorporated by 

reference in its entirety. 

This document is related to co-pending and commonly assigned U.S. 
15 Patent Application Serial No. 10/195,838, entitled "USE OF CURVATURE 

BASED FEATURES FOR BEAT DETECTION," inventor Sweeney, filed on 

July 15, 2002, the specification of which is hereby incorporated by reference in 

its entirety. 

This document is related to co-pending and commonly assigned U.S. 
20 Patent Application Serial No. 10/607,820, entitled "SIGNAL COMPRESSION 
BASED ON CURVATURE PARAMETERS," inventor Sweeney, filed on June 
27, 2003, the specification of which is hereby incorporated by reference in its 
entirety. 



25 Technical Field 

This document relates generally to cardiac rhythm management and in 
particular, but not by way of limitation, to detection and discrimination of 
arrhythmias in an implantable medical device. 



30 Background 

Implantable cardiac rhythm management devices typically monitor and 
process cardiac signals to provide therapy to the heart. Therapy may include 
delivering a pacing pulse to trigger a contraction of the heart or delivering a 
defibrillation shock to interrupt an abnormal heart rhythm. Some cardiac rhythm 
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* 

management devices also monitor cardiac performance or other physiological 
parameters for use in controlling the delivery of pacing pulses. 

Before delivering therapy, some implantable cardiac rhythm management 
devices are programmed to verify that an arrhythmia is occurring. In addition, 
5 therapy can be tailored to the patient's needs if the rhythm type is known. 

There exists an unmet need for providing improved arrhythmia detection 
and discrimination systems, devices, and methods to allow better diagnosis or 
treatment of patients. 

Summary 

1 0 The present subject matter provides a curvature based method of 

selecting features from a sampled signal. In one embodiment, the sampled 
signal includes a cardiac signal or an electrocardiogram. According to one 
aspect of the present subject matter, characteristic points are selected by sensing 
a cardiac signal on a real-time basis and upon receipt of each sample, computing 

1 5 curvatures on a continuous basis. Characteristic points are sometimes referred to 
as significant points or data points. 

In one embodiment, each "turn" in the signal is represented with a 
characteristic point that has measures of time, value, area and width. Time refers 
to a time denoting the center of the turn in the curvature series. Value refers to 

20 the amplitude of the input signal at the time of the center of the turn. Area refers 
to a value denoting the direction and degree of the turn. Width refers to a value 
denoting the duration over which the turn occurs. 

By the methods presented herein, each new characteristic point is known 
as soon as the turn which it describes has ended. In one embodiment, the stream 

25 of characteristic points are analyzed continuously as each new characteristic 

i 

point occurs. In one embodiment, the characteristic points are saved into a 
> buffer for analysis at a later time. 

All turns in an input signal are not the same. Some turns may be large 
deflections with sharp angles (for example, turns in a QRS signal) while other 
30 turns might be slight and not associated with big deflections (for example, noise 
in the signal). In one embodiment, the present subject matter includes a scheme 
to select prominent characteristic points in a signal and discard the smaller ones 
that may represent noise or only slight signal deflections. 



WO 2005/002669 PCT/US2004/020723 

The area for the characteristic point is an indirect measure of the angular 
turn in the signal. The maximum turn is approximately 180 degrees which 
represents a very rapid signal reversal. A large deflection in the signal, such as 
the R-wave (a full signal swing in 10 samples over, a period of 50 ms) would 
5 have large turn angle of 178 degrees. A small noise spike (such as a 10% full 
scale swing for a single sample) may also have a large turn angle of 175 degrees. 
The present subject matter treats these signal deflections differently. The R- 
wave has a large area whereas the noise spike would have a small area since it 
occurred over a lower number of samples. Thus, the characteristic point area is a 
1 0 surrogate for determining which characteristic points are important and which 
are not. 

The effect of the foregoing is that small characteristic points at or near 
the signal baseline value are ignored. According to one embodiment, an 
absolute area threshold value in the range of 0.5 is used for a rate-sensing 
15 channel. 

In one embodiment, the present subject matter automatically adjusts the 
threshold value to better adapt to varying patients and signals. 

In one embodiment, each characteristic point includes a time component, 
an amplitude component and an area component. According to the present 

20 subject matter, the time, amplitude and area components are derived from a 
curvature calculation performed using the sampled signal. 

Typical cardiac rhythm management devices determine the heart rate by 
sensing beats in a cardiac electrical signal and detecting when an R wave occurs 
in the signal. An R wave (sometimes referred to as a QRS complex) represents 

25 the major electrical event for each beat of the heart. In one embodiment, the 
plurality of characteristic points are analyzed to determine a heart rate without 
relying on beat sensing. In one embodiment, a first estimate of heart rate is 
generated by autocorrelating a particular function of characteristic points. In one 
embodiment, a time difference autocorrelation function yields the heart rate. In 

30 one embodiment, a second estimate of the heart rate (or beat frequency) is 

generated by a secondary process and the two estimates are reconciled to derive 
a single value for the heart rate. 



3 
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In one embodiment, an atrial series of characteristic points are generated 
for an atrial channel and a ventricular series of characteristic points are generated 
for a ventricular channel. The relative timing of the atrial series and the 
ventricular series are compared and the comparison forms the basis for 
5 discriminating among different rhythm types. Other methods are also 

contemplated including discriminating on the basis of variations in characteristic 
point timing and cross-correlations in the characteristic point domain. 

In one embodiment, a window is established before each ventricular 
characteristic point. Relative timing is evaluated by counting atrial and 
10 ventricular characteristic points relative to the windows and graphically plotting 
the results on a grid. The grid is divided into regions corresponding to different 
rhythm types. 

Other aspects of the invention will be apparent on reading the following 
detailed description of the invention and viewing the drawings that form a part 
1 5 thereof. 

k 

Brief Description of the Drawings 

In the drawings, like numerals describe substantially similar components 
throughout the several views. Like numerals having different letter suffixes 
20 represent different instances of substantially similar components. 

Fig. 1 illustrates an arbitrary signal as a function of time. 

Fig. 2 illustrates osculating circles fitted to an arbitrary signal as a 
function of time. 

Fig. 3 illustrates a signal expressed in analog-to-digital converter value as 
25 a function of sample number. 

Fig. 4 illustrates a signal expressed in analog-to-digital converter value as 
a function of sample number and having a constant curvature. 

Fig. 5 illustrates lobes of a curvature signal as a function of time. 

Fig. 6 graphically illustrates a method for determining area under a 
30 curvature lobe. 

Fig. 7 illustrates a characteristic point domain autocorrelation plot of 
characteristic points according to one embodiment of the present subject matter. 
Figs. 8 A, 8B, 8C, 8D, 8E and 8F illustrate a method to form a 

4 
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characteristic point time difference function according to one embodiment of the 
present subject matter and the resulting time domain autocorrelation. 

Fig. 9 illustrates a characteristic point domain autocorrelation plot 

* 

derived from a series of characteristic point. 

5 Fig. 10 illustrates fraction peaks in characteristic point domain from sub- 

harmonics. 

Fig. 1 1 illustrates a time function that represents intervals between 
characteristic points before and after a time shift. 

Fig. 12 illustrates an X-Y plane corresponding to a time relationship 
1 0 between atrial characteristic points and ventricular characteristic points for 
selected rhythm types. 

Fig. 13 illustrates a method of discriminating among rhythms according 
to one embodiment of the present subject matter. 

Fig. 14 illustrates a block diagram of an implantable device according to 
1 5 one embodiment of the present subj ect matter. 



Detailed Description 

In the following detailed description, reference is made to the 
accompanying drawings that form a part hereof, and in which is shown, by way 

20 of illustration, specific embodiments in which the invention may be practiced. 
These embodiments are described in sufficient detail to enable those skilled in 
the art to practice the invention, and it is to be understood that the embodiments 
may be combined, or that other embodiments may be utilized and that structural, 
logical and electrical changes may be made without departing from the scope of 

25 the present subject matter. The following detailed description is, therefore, not 
to be taken in a limiting sense, and the scope of the present invention is defined 
by the appended claims and their equivalents. 

. This document discloses, among other things, systems, method and 
apparatus involving cardiac rhythm management systems used for sensing a 

30 cardiac signal, including intracardiac electrogram or surface electrocardiogram 
(ECG). Such cardiac rhythm management systems include, but are not limited 
to, pacemakers, CRT devices, cardioverter/defibrillators, pacer/defibrillators, and 
drug delivery devices. However, it is to be understood that the present methods 



WO 2005/002669 PCT/US2004/020723 

and apparatuses of compressing a signal can be applied to heart beat detection as 
well as other signals related to cardiac activities, including, but not being limited 
to, mechanical motion, sound, pressure, acceleration, and impedance signals. 

5 Curvature Based Analysis 

For an arbitrary voltage versus time signal, the salient characteristics can 
be captured in terms of those points along the signal where it makes 
characteristic "turns." In general, these points are different from those that 
would be selected using first or second derivative criteria as shown in Fig. 1. In 

10 the figure, which illustrates arbitrary electrogram 100, the R wave, denoted as 
point 101, has a high second derivative dV/dt 2 . Further, the electrogram 
sections on either side of point 101 have high positive and negative slopes 
detectable with a first derivative criteria. However, those points with these large 
(or even maximum) slopes do not convey any particularly significant description 

15 of the electrogram. For example, at each time along a segment between points 
101 and 102, the waveform has a large negative slope but no point along this 
segment stands out significantly from any other. 

On the other hand, points 102, 103, 104 and 105 are neither maximums 
nor minimums of the electrogram or its derivatives. These points are descriptive 

20 of the arbitrary waveform shape. Points 102, 103, 104 and 105 are salient 

because they mark locations where the signal makes significant turns. The turn 
at point 101 is very sharp and the turns at points 103 and 105 are less sharp turns 
and more broad. The turns at points 102 and 104 are even less sharp but rather 
local. The present subject matter detects points 101, 102, 103, 104 and 105 

25 using a criteria based on the curvature of the signal. 

According to one embodiment of the present subject matter, signal 
curvature is illustrated in Fig. 2. In the figure, arbitrary electrogram 120 is 
shown. At each point along electrogram 120, an osculating circle can be found 
that fits tangentially with the local portion of the electrogram. Each circle has a 

30 radius. The curvature of electrogram 120 at that point is the inverse of that 
radius so that small circles have large curvatures while large circles have small 
curvature. Fig. 2 shows these circles at selected points. At point 122, the circle 



6 



WO 2005/002669 PCT7US2004/020723 

is larger, and thus the curvature is smaller, than at points 121, 123 and 124. At 
points 121, 123 and 124, the turns are sharper and the curvature is larger. 

In general, the curvature at point (X, Y) of an arbitrary curve in a two- 
dimensional space is expressed as: 

5 Curvature = d 2 Y/dX 2 

[l + {dY/dX} 2 ] 3 ' 2 

i 

As seen, the curvature is a non-linear combination of both the first and 
second derivatives of the curve. At those points along the curve where the first 

10 derivative is zero (for example, point 101 in Fig. 1), the curvature is equal to the 
second derivative and at points where the second derivative is zero (for example, 
any straight section regardless of its slope), the curvature is zero. 

The present subject matter calculates curvature of an arbitrary input 
signal on a sample-by-sample basis. 

15 Consider the question of dimensionality of the curvature for electrogram 

signals. When both X and Y have the same dimensions (say length) then 
curvature has a dimension of 1/length. For an arbitrary signal having a voltage 
as a function of time V(t), such as an electrocardiogram, the signal is 
transformed into a time-versus-time signal T(t) according to T(t) = V(t)/U where 

20 U is a constant with dimensions of voltage/time. With this transformation, the 
first and second derivatives of T(t) become 



25 



dT_(t) = dV_(t) • 1 / U which is dimensionless and 

dt dt 



d 2 T (t) = d 2 V ft) ♦ 1 / U which has dimensions of 1/time 
dt 2 dt 2 

30 and thus the curvature has dimensions of 1/time. Curvature is then expressed as: 

Curvature = d 2 TftVdt 2 - d 2 VftVdt 2 /U 

[1 + {dT(t)/dt} 2 ] 3/2 [1 + {dV(t)/dt/U} 2 ] 3/2 

35 which has dimensions of 1/time and U has a numerical value. 

Consider next, curvature as a function of signal gain or amplitude of the 
input signal. Assume arbitrary gain G is applied to the input signal to find a new 
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input signal F(t) wherein F(t) = G .T(t) = V(t) . G/U. The curvature of the 
gained signal is then: 

Curvature = d 2 F<Wdt 2 = d 2 Vffl/dt 2 . G/U . 

[1 + {dF(t)/dt } 2 ] 3/2 [1 + {dV(t)/dt . G/U} 2 ] 3/2 

5 

> 

The ratio G/U can be expressed as W having dimensions of time/voltage. 
In one embodiment, the input signal is a voltage sampled with a 12-bit analog- 
to-digital converter (ADC) having numerical voltage values in the range +/- 
2048 where each value represents a number of basic amplitude units AV, or 

1 0 voltage resolution. In one embodiment, the amplifiers in the present subject 
matter are adjusted so that the samples from V(t) largely fill this range. 
Furthermore, it is assumed that the samples are taken at a fixed rate and thus the 
time is represented by an integer number of samples with each sample 
representing a time interval AT = time resolution = l/(sample rate). 

1 5 Fig. 3 illustrates numerical X-Y space 130 which is denoted by sample 

numbers on the X-axis and the analog-to-digital conversion value on the Y-axis. 
Consider square 135 in X-Y space 130 which is 256 steps along the X-axis and 
256 steps along the Y-axis and rectangle 140 in voltage-time space 145 that this 
square represents. Rectangle 140 has a width of 256/(sample rate), in seconds 

20 and a height of 256/(voltage resolution), in volts. Depending on the values 
selected for AV and AT, this rectangle may or may not be a square in voltage- 
time space. 

In one embodiment, W is selected to require that a square in voltage-time 
space be represented by a square in sample-sample space. Under that condition, 

25 the curvature versus time relationships that exist in voltage-time space are 
preserved in sample-sample space. 

In Fig. 4, an artificial voltage-time curve is shown that forms a semicircle 
inside square box 150. Traversing from left to right across box 150 in voltage 
time space, the curvature of the signal is zero until the semicircle is encountered. 

30 The curvature jumps to a constant value, equal to 1/radius of the circle, and 

again jumps to zero at the end of the semicircle. For a particular value of W, the 
representation of this signal in sample-sample space 155 also has a constant 
curvature and for other values of W, the curvature would not be constant in 
sample-sample space. 

8 
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To maintain this relationship between voltage-time and sample-sample 
space, W is selected as follows. In voltage-time space, the box T wide by V high 
is taken to be square. In sample-sample space, the box is T/AT time-samples 
wide and GV/AV voltage-samples high and through a transformation, the voltage 
5 sample is converted into time-samples using U so that the sample-space square 
has dimensions of T/AT by VG/AV/U. For the box to be square in sample-space, 
and assuming the box (V by T) in voltage-time space is square, then (V/AV) . 
G/U = T/AT or G/U « W = (T/AT)/(V/AV). 

According to one embodiment of the present subject matter, curvature is 
10 based on the first and second derivative of the signal A least square cubic 

polynomial fit is used to reduce the noise that would otherwise result from using 
numerical estimates for the derivatives and using non-linear calculations to find 
curvature. 

The sampled voltage signal is expressed as V(t) = V(I • AT) where t=i • 
15 AT so that the fit of size N uses 2N+1 voltage samples centered on the time t, 
thus: 

V([i-n].AT), . . .V([i-2].AT), V([i-1].AT), V(i.AT), V([i+1].AT), V([i+2].AT), 
. . .V([i+n].AT). 

20 

This set of N sampled data points is used to make the least-squares cubic 
fit which is given as [i.AT + dt] = Ai + Bi . [dt] + Ci • [dt] 2 + Di . [dt] 3 where 
Ai, Bi, Ci, and Di are coefficients determined by minimizing the square error for 
the fit and dt represents a time step away from i*AT at which V^t is evaluated. 
25 The coefficients in the polynomial are denoted with an i to indicate that they are 
valid for the point i AT. 

Using the above equations, the curvature on a point to point basis, at the 
time i»AT, becomes 

* 

Curvature ( i.AT) = d 2 WtVdt 2 . W = 2 . Ci . W 
30 [l + {dV(t)/dt.W} 2 ] 3/2 [l + {Bi.W} 2 ] 3/2 

Sample points, however, will not necessarily fall at the times where the 
signal curvature has a maximum or minimum value. Thus, in one embodiment, 
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the curvature signal is integrated between adjacent sample points using estimates 
for the first and second derivative of the signal at the sample points. 

Further simplification yields an expression for average point curvature 

as: 

5 , Cava = 1 • {C^ 2 -Cf Bj»BM -q +C U + CgBfcuCfctfia^Bffq -C a } ^i^j., 

(BrB,,) 2 (B^l)™ (B M 2 +i) r/2 
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Cavgj = JL • c » + c *-' , Bj =Bj., 
2 (Bi 2 +1) 3/2 



As noted, curvature is computed on a sample-by-sample basis from the 
input signal. 

Consider next, a procedure for finding the characteristic points in the 

signal. 

1 5 Turns in the original signal are reflected as excursions above and below 

zero in the curvature signal. As suggested by curvature 160 of Fig. 5, each lobe 
above zero (for example, lobe 165) or below zero (for example, lobe 170) then 
represents a single turn in the input signal. Curvature lobes of opposite 
directions reflect opposite turns (leftward or rightward) in the signal. The area 

20 under each lobe reflects the total angle included in the turn. 

A point-by-point method is used to identify the lobes as they occur and to 

* 

find the area and centroid of each lobe. , 

Curvature values are generated on a sample by sample basis. At each 

sample time, CRV represents the current curvature value and CRV 0 id represents 
25 the previous value that is retained from the previous sample. When the 

curvature value is zero, the original signal is not turning so that no characteristic 

points may exist. The value CRV will rarely be exactly equal to zero. 

In one embodiment, a dead-zone surrounding zero is defined where the 

calculated curvature may be treated as equal to zero. For times where the 
30 curvature is within the zone, the signal is not turning significantly. The zone is 

defined with a curvature threshold value and extends above and below zero from 

+CRV t hresh tO -CRV t hresh. 

In one embodiment, the present subject matter identifies nine cases for 
considering the value of CRV relative to the threshold and the absence or 
35 direction of a lobe. These cases can be described as follows: 
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Case 1 : CRV>CRVthrEsh and not in a lobe. 

Here, the current curvature value is above the dead-zone and the 
curvature signal is not currently in a lobe. Thus, a positive lobe has just started 
so the positive lobe initialization calculations will follow. 

5 

Case 2: CRV th resh >CRV > -CRVthresh and not in a lobe. 

Here, the current curvature values is inside the dead-zone and the 

* 

curvature signal is not currently in a lobe. 

1 0 Case 3 : CRV<- CRV thresh and not in a lobe. 

Here, the current curvature value is bejow the dead-zone and the 
curvature signal is not currently in a lobe. Thus, a negative lobe has just started 
so the negative lobe initialization calculations will follow. 

1 5 Case 4: CRV>CRVthresh and in a positive lobe. 

Here, the current curvature value is above the dead-zone and the 
curvature signal is in a positive lobe. Thus, the positive lobe continuation 
calculations will follow. 



20 Case 5: CRV^n, >CRV>- CRV thre sh and in a positive lobe. 

Here, the current curvature value is inside the dead-zone and the 
curvature signal is in a positive lobe. Thus, a positive lobe has just ended so the 
positive lobe finalization calculations will follow. 



25 Case 6: CRV< -CRVthresh and in a positive lobe. 

Here, the current curvature value is below the dead-zone and the 
curvature signal is currently in a positive lobe. Thus, a positive lobe has just 
ended and a negative lobe has just started so that both the positive lobe 

4 

finalization calculations and the negative lobe initialization calculations will 
30 follow. 



Case 7: CRV >CRV,h TO h and in a negative lobe. 
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Here, the current curvature value is above the dead-zone and the 
curvature signal is in a negative lobe. Thus, a negative lobe has just ended and a 
positive lobe has just started so that both the negative lobe finalization 
calculations and the positive lobe initialization calculations will follow. 

5 

Case 8: CRVthresh >CRV> -CRV t hresh and in a negative lobe. 

Here, the current curvature value is inside the dead-zone and^e 
curvature signal is in a negative lobe. Thus, a lobe has just ended so the negative 
lobe finalization calculations will follow. 

10 

Case 9: CRV< -CRVthresh and in a negative lobe. 

Here, the current curvature value is below the dead-zone and the 
curvature signal is in a negative lobe. Thus, the negative lobe continuation 
calculations will follow. 

15 In one embodiment, hysteresis is used to aid in identifying curvature 

lobes. After having started a lobe, one embodiment provides that the curvature 
value cross a threshold value closer to zero in order for the lobe to finish. Thus, 
hysteresis introduces another threshold value. 

Consider next the metrics for characterizing each lobe according to the 

20 present subject matter. In one embodiment, those metrics include the total area 
of the lobe, the time of the lobe area centroid and the value of the original data at 
the time of the area centroid. In one embodiment, other metrics are used, 
including, for example, the peak curvature in the lobe, the time of the peak 
curvature, and the times of the lobe start and lobe finish. 

25 Fig. 6 shows an example of curvature lobe 185 as a series of calculated 

curvatures at each sample time. Tick marks 1 80 along the X axis represent the 
actual samples and thus are separated by AT. Zero curvature is shown as the X 
axis. The threshold curvature values, +CRV th resh 190 and -CRVthresh 195 are 
shown as horizontal lines above and below the X-axis. Curvature points 200- 

30 245 are shown as small solid circles. 

As shown in the figure, initial curvature values at curvature points 200 
and 205 are within the dead-zone of the threshold curvature values and thus, no 
lobe is yet established. Curvature point 210 is above the threshold and 

12 
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corresponds to Case 1 described above. The magnitude of initial area 250 is 
calculated as soon as the lobe is started according to case 1 described above. 
Upon determining curvature point 215, a contribution to the magnitude of main 
area 255 is calculated according to case 4 described above. In addition, the 
5 magnitude of main area 255 is increased with the determination of curvature 
points 220, 225, and 230. Upon deteraiining curvature point 235 lying below 
+CRVthresh 190, the magnitude of final area 260 is calculated for the region 
below curvature lobe 1 85 and between curvature point 230 and the intersection 
with +CRVth res h 190. Curvature points 235, 240 and 245 lie in the dead-zone and 

■ 
■ 

10 do not contribute to an area calculation. The area of curvature lobe 185 is 
calculated to include the sum of initial area 250, main area 255 and final area 
260. 

In one embodiment, a value M of a counter is maintained by a processor 
to monitor sample-by-sample development of a lobe. When the 'start lobe* 

1 5 calculations are performed, the value for M is reset to zero and incremented with 
each subsequent point that remains outside the dead-zone. Values for M are 
illustrated at the different curvature points in Fig! 6. 

To find the area under the curvature lobe, the value of M is set to zero 
and the initial area of a curvature lobe is computed as initial area = (Q + 

20 CRVth re sh) • (Ci - CRVttosh) / (CrCi-i) when the lobe is started. For each 

successive curvature value, the value of M is incremented and the additional area 
contribution is computed as additional area = - (M-l) • Q + M • Cm. When the 
lobe ends, the value of M is not incremented and the final area contribution is ' 
computed as final area = M • C M + (Q + CRV threS h) • (Q - CRVthreshVCQ-Ci-i). 

25 The sum of the initial, final, and all additional areas may then be multiplied by 
AT/2 to find the area under the curvature lobe. 

■ 

To find the first moment of the area under the curvature lobe, a similar 
approach is used. When the initial area of the curvature lobe is computed, the 
initial moment of the area is also computed as initial moment = - {Ci + 
30 2 CRV lhr esh} • Y 2 where y = (CRV 

thresh " Q-i)/(Cj-Ci.i). For each successive 
curvature value in the lobe, an additional moment contribution is computed 
using the same M values as for the area computations. This additional moment 
contribution is computed as additional moment = {Q • (3M-1) + Cm • (3M-2)}. 
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When the lobe ends, a final moment contribution is computed at the same time 
that the final area is computed. The final moment is computed as final moment 
=[CRVth reS h • { 3M + 27} + Cm • {3M + 7} ] • 7. 

The sum of the initial, final, and all additional moments may then be 
5 multiplied by AT 16 to find the first moment of the area under the curvature lobe. 
The time of the centroid of the curvature lobe area is found by dividing the first 
moment of the curvature lobe area by the curvature lobe area as time = first 
moment of area/area. This time is with respect to the time of the curvature point 
that started the lobe. 

1 0 For both the area and the first moment of the area, the calculations are 

constructed in such a way that the contribution of each curvature sample may be 
summed into a total area or total first moment of the area as the samples are 
collected. Thus, for all calculations, the factors include the present and previous 
curvature values, the threshold values, and the counter M. In this manner the 

15 area and the time of the centroid of the area are generated as soon as the lobe is 
ended. 

In one embodiment, the present subject matter identifies characteristic 
points generated from point-by point processing of an input signal, hi one 
embodiment, the characteristic points in the signal are detected and saved into a 

20 buffer. In one embodiment, each characteristic point is a set of values including 
the time of the characteristic point, the value of the input signal at that time, and 
a value describing the direction and extent of curve in the signal that produced 
the characteristic point. 

Consider next a method of extracting characteristic points. In one 

25 embodiment, characteristic points of an electrocardiogram are extracted using a 
sampling at 400 Hz. In one embodiment, sampling is performed at 200 Hz. The 
sampled data is filtered using a 5-point (400 Hz) or 3-point (200 Hz) running 

■ 

average filter and a five least squares regression average-point curvature method 
to select characteristic points relating to the QRS complex. 
30 For detecting curvatures associated with slower morphologies, sampling 

is at a slower rate, such as, for example, 50 Hz. An example of a wave having 
slower morphologies would be the T-wave in an electrocardiogram. 
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In one embodiment, multiple sampling rates are used for characteristic 
point extraction. 

In one embodiment, dual-rate sampling is performed with sampling at 
200Hz providing fast characteristic points and sampling at 50 Hz providing slow 
5 characteristic points and a 3-point running average filter is used. Consider an 
example of dual-rate sampling using 200 Hz, yielding fast characteristic points, 
and 50 Hz, yielding slow characteristic points, with 5-point least squares 
regression size and 3-point running average filtering at both rates. 

In one embodiment, curvature is found by using 5-point least squares 
10 regression filtering. The coefficients for finding the linear (Bi) and parabolic 
(Ci) fit coefficients in the cubic least squares regression to the data (D^ 2 to Dj +2 ) 
at the point i are as follows: 



15 





Linear 




Parabolic 


P-2= 


(1/12) * Rate 


Q-2 


= (10/70)* Rate 2 


P-1 = 


(-8/12) * Rate 


Q-i 


= (-5/70)* Rate 2 


Po = 


0 


Qo 


= (-10/70)* Rate 2 
= (-5/70) * Rate 2 
= (10/70)* Rate 2 


Pl = 


(8/12) * Rate 


Qi 


P 2 = 


(-1/12) * Rate 


Q2 



20 where Rate is the sampling rate in sample/sec, Bi is the sum of the products of 

these P least squares regression coefficients multiplied by the 5 corresponding 

data points centered on i (Dj_ 2 to D i+2 ) and Ci is the sum using the Q least squares 

regression coefficients. Curvature is then computed as: 

Curvature = 2 • Ci * W 
25 [l + {Bi.W} 2 ] 3/2 

where W is a constant. Since W appears in conjunction with fit coefficients Ci 

or Bi, the value can be incorporated into the computations as a data signal gain 

such that D f (i) = W • D(i). Then, the curvature expression becomes 

30 Curvature = 2«C'i 

[l + {B'i} 2 ] 3/2 

where B'i and Ci are found using the D' data points rather than the D data points 

with the least squares regression coefficients. 

35 The expression for curvature may then be rewritten as: 

Curv' = Curvature . Gn = C'i»2.Gn . = C"i 

[l + {B'i} 2 ] 3/2 [l + {B'i} 2 ] 3/2 
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where Gn is an arbitrary value which is incorporated (with 2) into C". If the 
curvature thresholds used for detection of characteristic points is adjusted 
accordingly, it does not matter whether point-by-point curvature or Gn times the 
point-by-point curvature is calculated. Incorporating the values 2 . Gn into C"i 
5 can be accomplished by changing the least squares regression coefficients Q_2 
through Q2 into new least squares regression coefficients (Q 1 ^ through Q'2 ) by 
multiplication with 2Gn. If Gn is selected as 7/(6 Rate) and the term Fn = 
Rate/ 12 is incorporated as additional data signal gain, then the least squares 
regression coefficients become: 

10 Linear Parabolic 

P M - 2 = 1 Q n - 2 = 4 

P"-i=-8 5 Q"-i= -2 

P"o = 0 Q%= -4 

P f, i = 8 Q") = -2 

15 P" 2 = -1 Q» 2 = 4 

with D"(i) = 2 . W . Fn . D( i ); and with Ci" = Q"_ 2 • D M (i-2) + Q".,-. D"(i-1) + 

Q"o • D"(i) + Q\, . D"(i+1) + Q" +2 • D"(i+2); and with Bi" = P"_ 2 • D"(i-2) + P". 

1 .D"(i-l) + P ,, o-D ,, (i) + P"-i •D , '(i+l) + P" +2 .D"(i+2). 

Note that the coefficients are all powers of two (either 0, 1, 2, 4 or 8), 

20 thus simplifying hardware or firmware implementation of the present subject 

matter. 

When using two sample rates for determining characteristic points, the 

computed curvatures will be different when computed at different rates. Thus, 

two threshold levels and parallel sets of computations are used. Also, the 
25 running average filtering is different with the fast rate filtering using five 

consecutive points while the slow rate filtering uses every fourth point spaced 

over sixteen points. 

In one embodiment, the linearity of the least squares regression operation 

is used. The running average filtering operation are performed at the fit 
30 coefficient level (i.e. the Bs and Cs) rather than the data level. In this 

embodiment, two filtering operations are performed (one for Bs, one for Cs) at 

each rate. 

In one embodiment, the fast and the slow running average filtering 
operations are performed in parallel. 
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In one embodiment, the following method is used to select the values for 
the characteristic points. 

Starting at the time of one sample, the sensing hardware and software of 
the present subject matter begins acquiring the next sample so that the actual 
5 sample reflects the nature of the signal between the last sample and the current 
sample. At the actual time of the sample, the value of the sample may be 
ascribed to the time half-way between the current sample and the previous 
sample. Samples at the slow sample rate apply to the signal two fast sample 
steps earlier than the actual sample time. 

10 In one embodiment, a circular data buffer is used. Data is acquired at the 

fast sample rate and used to fill a circular data buffer. Once the timing of a fast 
or slow characteristic point is determined, the data value for that characteristic 
point is found from the data stored in the circular buffer. Thus, the data buffer is 
sized so that the requisite data is still available in the buffer. In one embodiment, 

1 5 the buffer is sized to hold about 500 millisecond of data or 256 data points. In 
one embodiment, the circular buffer is implemented as a memory array with an 
index pointer. In one embodiment, the index pointer includes a binary counter 
and the buffer size is a power of two. 

In one embodiment, running average filtering (RAF) is performed on the 

20 input data. In various embodiments, the filtering is provided by a circular buffer 
or a shift register set. The values in the circular data buffer reflect the fast 
running average filtering. In one embodiment, a hybrid approach is used where 
the fast running average filtering is performed and then the data is placed into a t 
circular buffer. 

25 In one embodiment, running average filtering reduces curvature noise 

relative to the signal content of fast curvature signals. An electrocardiogram 
signal typically takes relatively large excursions away from zero curvature at the 
times of fast turns in the signal. However, when the fast turns are less severe, 
the curvature signal does not move as far away from zero and curvature noise 

30 can prematurely end the curve lobe. The result is that one small curvature lobe 
is broken into two or more even smaller lobes. Thus, in one embodiment, the 
curvature noise relative to the signal content is reduced, for fast curvature 
signals, by running average filtering. 
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Filtering of the input signal does not achieve the same effect as filtering 
the actual curvature signals. In one embodiment, 3-point running average 
filtering is applied to the fast curvature values. 

In one embodiment, fast curvature values are computed as average point 
5 curvature and 3-point running average filtering of the curvatures so that the fast 
curvature value computed when the current sample is taken applies to a time in 
the input signal that is five time steps earlier than the current sample. 

In one embodiment, slow curvature values are computed at 1/4 the fast 
sample rate and are computed as average-point curvatures so that when the every 
10 fourth fast sample is taken, slow curvature is computed and it applies at a time in 
the input signal that is four slow (sixteen fast) time steps earlier than the current 
sample. 

In one embodiment, the detection of a curvature lobe and characteristic 
points entails both fast and slow characteristic point detection performed in 
15 parallel. In one embodiment, separate curvature thresholds and curvature area 
thresholds are used for the fast and slow operations. 

In one embodiment, a curvature lobe is ended when the curvature signal 
crosses the same threshold used to start the lobe but in the opposite direction of 
the crossing that started the lobe. In one embodiment, a hysteresis is used so that 
20 a curvature lobe is ended when the curvature signal crosses a threshold that is 
smaller (i.e. closer to zero) than the threshold used to start the lobe but in the 
opposite direction of the crowing that started the lobe. In one embodiment, a 
curvature lobe is ended when the curvature signal crosses zero as the hysteresis 
value. 

25 In one embodiment, the value of W is selected according to the following 

procedure. Select an approximate range for W from the following analysis. At a 
200 Hz sampling rate, the QRS takes approximately five samples. The signal at 
the peak is approximated by 128. If the QRS were to be approximated (in 
curvature space) as a semicircle arching from zero to the peak then back to zero 

30 in five points, then the radius of that semicircle would be 0.01 sec and the peak 
would have a value of 0.01. Thus, the value of W would be the one that, when 
multiplied, the peak (128) would have the value of 0.01, i.e. W = 0.01/128 = 
0.000078. In one embodiment, W = 0.000122= 1/8192 which is a power of two. 
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Consider selection of curvature threshold values for characteristic point 
detection. Thresholds, in one embodiment, reduce or eliminate the noise 
characteristic points and computational burden. This noise is the natural results 
of the curvature signal hovering around zero (or some other value) when the raw 
5 signal does not otherwise have characteristic turnings. 

In one embodiment, for a curvature lobe to be detected, the point-wise 
curvature falls outside of a dead-band formed by the threshold values on either 
side of zero. Reducing this threshold towards zero increases the number of 
curvature lobes that are detected and increases the computational burden. In one 

10 embodiment, a threshold value is selected by using approximately a 1 minute 
epoch of the input signal to make a histogram of the fast curvature values 
computed at each sample time. The threshold value is selected as one which 
corresponds to 20% of the peak. This value is then used as the threshold for both 
the fast curvature and slow curvature feature selection process. 

15 In one embodiment, the noise lobes are removed by requiring that a 

detected lobe must exceed a critical area limit. The noise lobes have generally 
small areas while the real lobes associated with characteristic turns in the signal 
usually have substantial areas. 

In one embodiment, different values are used for the fast curvature and 

20 slow curvature area thresholds. Using the number of beats occurring in the one 
minute epoch used to set the gain and curvature thresholds, select a desired 
number of fast and slow characteristic points to find for these beats. For 
example, in one embodiment, five fast characteristic points per beat and eight 
slow characteristic point per beat are used as targets. The area thresholds are 

25 then adjusted until these target numbers of characteristic points are found. In 
one embodiment, a lower area limit of 0.1 is used. 

In one embodiment, a template is established for classifying beats. For 
example, a characteristic point template includes a rectangular box covering a 
time (relative to a fiducial characteristic point) and signal amplitude. The 

30 template also includes a curvature area sign. 

Rate Estimation 

Heart rate, expressed in units of beats per minute, can be determined 
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without relying on detection of individual beats. Mathematically, the heart rate 
is determined by the number of characteristic points per minute divided by the 
number characteristic points per beat. The calculation of the number of 
characteristic points per minute is straightforward. An autocorrelation function 
5 is evaluated to determine the number characteristic points per beat. 

In one embodiment, the number of characteristic points per beat is 
calculated by performing the autocorrelation function in the time domain. In one 
embodiment, the number of characteristic points per beat is calculated by 
performing the autocorrelation function in the characteristic point domain. 

10 In the series, or stream, of characteristic points describing the input 

signal, let the individual characteristic points be denoted as CPi, CP 2 , CP 3 , 
CP n . Each characteristic point has a time value, an area (turn area) value and a 
voltage value, and for the i* characteristic point, these are denoted as CPi T> CPi A 
and CP,- v, respectively. 

1 5 The values for the autocorrelations range between +1 and -1 . The shift 

index (k) for the autocorrelations denotes the shift in the characteristic point 
domain. For example, at k=0, there is no shift and at k=l, the characteristic 
points are shifted by one. 

In one embodiment, the autocorrelation of characteristic points in the 

20 characteristic point domain entails, for each value of the shift index k, 

multiplying like metrics of corresponding characteristic points and summing to 
produce a value. For example, when k=3, the value of the series includes the 
sum of the product of CPiv and CP 4V , CP 2V and CPsv, CP 3V and CP 6V and so on. 
The results are graphically depicted in the autocorrelation plot of Fig. 7. 

25 In the figure, the vertical axis is an autocorrelation coefficient based on an 

autocovariance function divided by a variance function and the horizontal axis is 
the shift index k. The plot in the figure has a peak at a shift index of five which 
indicates that on average, each beat is represented with five characteristic points. 
Using exemplary data collected over an epoch having a duration of 6.288 

30 seconds, a total of 80 characteristic points were generated for a rate of 12.722 
characteristic points per second. With five characteristic points per beat, the 
heart rate is calculated as 152.67 beats per minute. 
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In one embodiment, an autocorrelation is performed using the time 
differences between successive characteristic points. For example the difference 
in time for adjacent characteristic points can be expressed as CPj dt - CPi t - 
CPi-i t. In addition, the difference in area for adjacent characteristic points can 
5 be expressed as CPi da = CP,- a - CPm a and the difference in voltage can be 
expressed as CPj dv = CPj v - CPm v. 

In various embodiments, the difference functions CPi dt, CPi da and CPi 
dv are autocorrelated to determine the number of characteristic points per beat. 
In various embodiments, CPiA, CP iV , CPj D T, CP iD A> and CPiov are autocorrelated, 
1 0 either alone or in combination (for example, the product CPiv x CPiDT x CP^), to 
determine the number of characteristic points per beat. 

By way of example, an autocorrelation of the characteristic point voltage 
values in the characteristic point domain can be expressed as: 

15 i = N i = N 

A v (k) = . E(CPi v - CPiv Av G )(CP i+k v - CP iV avg) / 2(CPj v - CP iv avg) 2 
1=1 i= 1 

where N is the number of characteristic points, CPiv avg is the average 
20 value of the N characteristic point values and k is the index in the characteristic 
point domain. ^ 

■ 

In one embodiment, heart rate is found from the autocorrelation of 
characteristic points in the time domain. In one embodiment, the cardiac signal 
is reconstructed for autocorrelation from the series of characteristic points. In 

25 one embodiment, the reconstruction produces evenly spaced samples for the 
autocorrelation. In one embodiment, the reconstruction and autocorrelation are 
performed in closed form using the characteristic point voltages and times from 
the series of characteristic points. 

In one embodiment, only timing information for the characteristic points 

30 in the cardiac signal are used for the autocorrelation. In one embodiment, a time 
difference between successive characteristic point times in the series of 
characteristic points versus the time of each characteristic point in the series 
forms a time difference function. In one embodiment, the time difference 
function is constructed to have evenly spaced samples for autocorrelation. In 

35 one embodiment, the autocorrelation of the time difference function is performed 
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in closed form using only the characteristic point times from the series of 
characteristic points. In general, the first peak of the autocorrelation is equal to 
the reciprocal of the heart rate. 

Fig. 8 A illustrates a selected epoch of cardiac signal 300 having 

5 characteristic points 305 marked thereon. The epoch includes five heart beats, 
each having between four and six characteristic points. Fig. 8B illustrates the 
amplitude of the signal, at each characteristic point, as a function of time. Fig. 
8C and 8D illustrate the time differences between successive characteristic 
points in the series of characteristic points. When plotted against the times of the 

1 0 characteristic points, these form a time difference function which may then be 

i 

autocorrelated. In Fig. 8C, the time difference 310 between adjacent 
characteristic points 305 A and 305B is transformed to an amplitude denoted by 
the horizontal line 315. Thus, the time difference function between points 305A 
and 305B is equal to this line segment 315. Fig. 8D illustrates the time 
1 5 differences for all characteristic points in the epoch. Fig. 8E illustrates the time 
difference after having removed the characteristic points. The waveform shown 
in Fig. 8E is depicted as an autocorrelation plot in Fig. 8F which shows that the 
first peak occurs at time 0.360 seconds which corresponds to a beat rate of 166 
bpm. 

20 In one embodiment, heart rate is found from the estimated number of 

■ 

characteristic points per beat and the times of the characteristic points in the 
series of characteristic points. In one embodiment, an average interval 
representing one beat is found and the interval is the reciprocal of the heart rate 
for the rhythm. For example, if the estimate for the number of characteristic 
25 points per beat is 5, then the average time interval between each pair of 
characteristic points that are 5 characteristic points apart is used to find an 
average time interval for the rhythm as follows: 

i = N 

30 Average Interval = 2 (CPj t - CPj- Kp t) 

i-Kp 

where Kp is the estimated number of characteristic points per beat. 
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Other autocorrelations are also contemplated. For example, in various 
embodiments, the characteristic point voltage values and area values are 
autocorrelated. 

Sub-harmonics in the autocorrelation function can be used to determine 

» 

5 the heart rate. The sub-hannonic frequencies are multiples of the fundamental 
frequency of the heart rate. For example, if the underlying rhythm is a 
monomorphic ventricular tachycardia (MVT), then the same basic beat shape is 

f 

repeated at regular intervals. In the characteristic point domain, this means that 
the autocorrelation plot will have multiple peaks. For example, an 

10 autocorrelation plot for a MVT rhythm exhibits a peak at a shift index of 5 as 
well as the sub-harmonics corresponding to the shift index at 10, 15, and 20. 

In one embodiment, the sub-harmonics are used to find the rate. In plot 
330 of Fig. 9, exemplary autocorrelations for DT are illustrated. The 
autocorrelations in the figure are derived from a particular cardiac epoch. In plot 

15 330, it is not clear if first peak 320 for the corresponding atrial cardiac epoch is 
at k=5 or k=6. The location of second peak 325, derived from the first sub- 
harmonic, can resolve the ambiguity. 

In one embodiment, a method to process the autocorrelation provides an 
improved resolution to find the peaks by incorporating sub-harmonics. As in 

20 Fig. 9, the autocorrelation has values at k=l, k=2, etc. Initially, additional points 
are added to the autocorrelation by adding points at 1/4 step resolution. For 
instance, in Fig. 9, additional points are added at 2.25, 2.5, 2.75, 3.25, 3.5, 3.75, 
and so forth. The autocorrelation values for these additional points are found by 
interpolating between the original points. Next, a sub-harmonic autocorrelation 

25 is found by averaging the values at harmonics. For example, at k=2, the value 
for the sub-harmonic autocorrelation is found by averaging the points at 2, 4, 6 
etc. of the autocorrelation. At the point k=2.25, the value for the sub-harmonic 
autocorrelation is found by averaging the points at k=2.25, 4.5, 6.75, etc. of the 
autocorrelation. Fig. 10 illustrates an example of the sub-harmonic 

30 autocorrelation from Fig. 9. In Fig. 10, the peak is now more clearly at k=5.5 
which is between the steps k=5 and k=6 of the original correlation. In one 
embodiment, the averages of the first three sub-harmonics are used in to form 
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the sub-harmonic autocorrelation. In one embodiment, the peak of the sub- 
harmonic autocorrelation is used to estimate the rate of the rhythm. 

. Autocorrelation in the time domain can be used for identifying 
periodicity in a signal. Periodicity can also be identified by performing the time 
5 domain autocorrelation of the voltage signal using the characteristic points 
voltage values. In one embodiment, the time domain signal is represented by 
straight line segments connecting consecutive characteristic points. The line 
segment representation of the signal is then autocorrelated using integration of 
line segments. The endpoints of the line segments are the characteristic points, 

10 and thus, the calculations are simplified. 

For example, consider a signal described by N characteristic points, 
where CPi is the oldest characteristic point, CP 2 is the next oldest and so on. A 
closed form autocorrelation for a given time-lag Q is then computed by sliding 
each line segment by Q and then integrating (for that line segment) the product 

1 5 of that line segment and the portion of the signal with which it overlaps. The 
values are obtained for all line segments to find the autocorrelation for the time- 
lag Q. The same operation is then performed for other values of time-lag to find 
the time domain autocorrelation (autocorrelation vs. time-lag). 

The heart rate can also be determined by autocorrelating a function of the 

20 characteristic points in the time domain. For example, in one embodiment, a 
closed form autocorrelation is calculated, as shown in Fig. 1 1 . 

Consider a sequence of characteristic point time differences CPj,T where 
there are N characteristic points. Time difference function 41 5 A is established 
and has a value, at each characteristic point, equal to the time interval between 

25 CPj and CPj+i as shown with the heavy black line and solid black circles in chart 
405. Time function 415 A is shifted in time by Q as represented as 41 5B in chart 
410. Each point along the shifted function 41 5B now falls onto some point of 
the un-shifted function 41 5 A. For example, point 420B (derived from shifted 
point 420 A) on chart 410, is projected onto point 420C on chart 405. Similarly, 

30 the values of the un-shifted function projects onto the time-shifted function. For 
example, point 425A of chart 405 projects to point 425B of chart 410. The 
autocorrelation of this function in the time domain is found by determining the 
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area under the product of the curves which is equivalent to adding the area 
contribution for each portion between the line vertical lines shown in Fig. 11. 

Next, consider the shifted curve 41 5B on chart 410. Each segment 
between the vertical lines is bounded in time by either a shifted or un-shifted 
5 characteristic point. Furthfer, these bounding endpoints have a value equal to 
either the shifted value or the un-shifted value respectively. Each segment also 
aligns with a second set of values (shown in chart 405) as either un-shifted (solid 
circles) or shifted (open circles) values. The values for the open circles (chart 
405) or open circles (chart 410) are found by shifting the characteristic points 

10 either forward (chart 405) or backwards (chart 410) in time. In one embodiment, 
the calculation of area is found from the shifted line segments, and the 
corresponding un-shifted line segments. For example, the line segment between 
420B and the next point on the shifted function 41 5B falls entirely on the same 
line segment after point 430 of the un-shifted function 41 5A. The values for 

15 420Band430 are known. Because 415B is a time difference function, the 
duration after 420B until the next point is also the value of 420B. Thus, the 
contribution to the autocorrelation integral for this line segment is equal to the 
value of point 420B squared times the value of the point 430. The 
autocorrelation is the sum of all line segments. 

20 The time domain autocorrelation of the characteristic points time 

difference function can be used for estimating rate. Various methods are used to 
estimate the rate. In some methods, the autocorrelation between 0.15 sec and 0.6 
seconds time lag is fit with a straight line via least squares regression and this 
trend-line is removed from the autocorrelation data in this same range. In one 

25 method, a point-by-point scan of the de-trended autocorrelation is performed to 
find all local peaks. Each local peak, is only then considered if it is faster than 
the first sub-harmonic of the first such peak. For example, if the first peak in the 
autocorrelation represents a rate of 300 bpm, then peaks less than or equal to 150 
bpm are not considered. All of the peaks remaining in consideration have a 

30 specific lag time and an autocorrelation value. In one embodiment, the 

autocorrelation values are used to weight the rates associated with those peaks to 
find an average rate for the rhythm. 
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In one method, the value of the de-trended autocorrelation function in the 
range from 0.15 to 0.6 seconds time lag and above zero are used to find an 
average time lag as SUM (corr*lag) / sum(corr). The average time lag is used to 
find the rate. In this method, the product of each de-trended autocorrelation 
5 value that is between 0.15 seconds and 0.6 seconds and greater than zero is 
multiplied by the corresponding time-lag and summed to find SUM(corr*lag). 
Also, the de-trended autocorrelation values that meet these criteria are 
themselves summed to find SUM(corr). 

In one embodiment, the average autocorrelation value associated with the 
10 peak is used to find the rate. 

Other methods are also contemplated for finding the heart rate. 

Rhythm Discrimination 

In one embodiment, analysis of characteristic points is used to 

1 5 discriminate among different rhythm types. 

Rhythm types are associated with a relationship between an atrial cardiac 
signal and a ventricular cardiac signal. Thus, to discriminate among different 
types, an atrial input signal and a ventricular input signal are used to generate a 
series of atrial characteristic points and a series of ventricular characteristic 

20 points. In various embodiments, the rhythm can be discriminated by estimating 
the variation in characteristic point timing and cross-correlations in the 
characteristic domain. In one embodiment, the relative timing of characteristic 
points in the atrial and ventricular rate channels is used to discriminate. 

In one embodiment, the individual characteristic points in the atrial 

25 channel are compared with those in the ventricular channel. 

In one embodiment, a time window is established in advance of each 
ventricular characteristic point. For an appropriately sized window, a 
characteristic point in the ventricular complex should have one or more atrial 
characteristic points within the window. Each ventricular characteristic point is 

30 considered individually and the number of atrial characteristic points that are 
contained in such a window is determined. In one embodiment, a constraint 
requires that each atrial characteristic point is counted only once. This count is 
referred to as a forward count and is denoted A: V and is equal to the number of 
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atrial characteristic points which fall in a time window of one or more 
ventricular characteristic point. 

The dimensions of the time window are selected to yield accurate rhythm 
discrimination results. In one embodiment, the time window is located between 
5 300 milliseconds and 100 milliseconds before a ventricular characteristic point. 
In one embodiment, the time window is located between 250 milliseconds and 
150 milliseconds before a ventricular characteristic point. Other dimensions are 
also contemplated. In one embodiment, the time window dimensions are 
selectable based on manual inputs or automatically based on a measured or 
10 calculated parameter. In one embodiment, a window dimension is remotely 
selected and communicated to an implantable device. 

Consider the resulting count, for a simple conducted rhythm where each 
atrial depolarization corresponds with a ventricular depolarization after some 
delay. In such a case, most of the atrial characteristic points will be counted so 
15 that the ratio of counted atrial characteristic points to the total number of atrial 
characteristic points (denoted as Ac nt ) will be nearly unity. 

In one embodiment, far-field sensing of ventricular depolarizations on 
the atrial channel may cause additional characteristic points in the atrial channel 
but these additional characteristic points would generally be simultaneous or 
20 later than the ventricular characteristic points so would not fall inside the 
windows. 

If the heart is beating in atrial fibrillation (A-fib), rather than a one-to-one 
conducted rhythm, then the total atrial characteristic points count will be larger. 
Also, a smaller percentage of atrial characteristic points, will fall into a 
25 ventricular characteristic point window since not all atrial complexes would be 
conducted and lead to a ventricular complex. Thus, the ratio should become 
smaller. 

If the heart is beating in MVT, polymorphic ventricular tachycardia 
(PVT) or ventricular fibrillation (VF), then the total number of atrial 
30 characteristic points would not increase but a smaller number of characteristic 
points would be associated with a ventricular characteristic point time window 
because of the loss of A-V synchrony. On the other hand, more ventricular 
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characteristic points would be present to increase the A: V number. Accordingly, 
the ratio is unlikely to be greatly affected. 

Next, consider the same timings from the point of view of the ventricular 
characteristic points. Each ventricular characteristic point is analyzed and a 
5 count is made of the number of ventricular characteristic points that have at least 
one atrial characteristic point within a time window for the ventricular 
characteristic point. This count is referred to as a backward count and is 
denoted V: A and is equal to the number of ventricular characteristic points 
which have one or more atrial. characteristic points in the advance time window. 
10 For a one-to-one conducted rhythm, most of the ventricular characteristic 

points will meet the criteria so that the ratio of counted ventricular characteristic 
points to the total number of ventricular characteristic points (V cn t) will be nearly 
unity. 

If the heart rhythm was actually A-fib instead of a one-to-one conducted 
1 5 rhythm, then the total number of ventricular characteristic points would be 
greater reflecting an elevated rate. However, since each ventricular complex 
would be from a conducted atrial complex, the percentage of ventricular 
characteristic points that are associated with atrial characteristic points would 
remain about the same even though there are more ventricular characteristic 
20 points. 

If the actual rhythm was an MVT or PVT/VF, then the total number of 
ventricular characteristic points would increase and a smaller number would be 
associated with an atrial characteristic point in their time window. Accordingly, 
the ratio would also be reduced. 

25 The ventricular and atrial characteristic points are found over the selected 

epoch. If the rhythm is one-to-one conduction, then there is the same number of 
beats in each channel. Without knowing the number of beats in the epoch, an 
estimate can be made for the number of atrial characteristic points per beat 
relative to the number of ventricular characteristic points per beat, denoted here 

30 asAcnt/Vcnt. 

If the rhythm is one-to-one, then the percentage of atrial characteristic 
points associated with ventricular characteristic points should approximate the 
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percentage of ventricular characteristic points associated with atrial 
characteristic points, however, the number of characteristic points per beat may 
be different in the two channels. 

In various embodiments, the value V:A/V cn t or the value [(A:V/Acnt) / ( 
5 V: A /V cn t) ] / [A^nt / Vent], which can be rewritten as (V cn t/Acnt) 2 /(A: V * V/A), is 
used to discriminate among rhythms. 

In one embodiment, a plot is generated to graphically illustrate the 
relationship between atrial characteristic points and ventricular characteristic 
points. In Fig. 12, a plot is shown for approximately 500 rhythms that had either 
10 supra-ventricular origin (SVT) or ventricular origin (VT). The X, or horizontal 
axis is 100*(V:A/Vcnt) and the Y, or vertical axis is 100*(A:V * Vcnt/Acnt 2 ). 

The different regions of the X-Y plane correspond to the different rhythm 
types. For example, the SVT rhythms fall generally in the range below about 
150 on the Y-axis and between 50 to 100 on the X-axis. In contrast, VT rhythms 
1 5 generally fall much closer to zero on the X-axis reflecting the expected strong 

■ 

decrease of V:A/V cn t. Furthermore, the VT rhythms fall higher on the Y-axis 
since V cm is generally higher. 

In one embodiment, an X and Y value for a rhythm is found by 
processing the timings of the atrial and ventricular characteristic points 

20 according to the present subject matter. In one embodiment, the position of the 
X and Y values on the plot is used for classifying the rhythm. In one 
embodiment, the classification is between either SVT or VT. In one 
embodiment, the classification is made by comparing the X-Y point with fixed 
boundaries on the plot. In one embodiment, the boundaries are comprised of one 

25 or more straight line segments. In one embodiment, the rhythm is classified as 
SVT if the X-Y point falls on one side of the boundary and is classified as VT if 
it falls on the opposite side of the boundary. 

Fig. 13 illustrates exemplary method 450 according to one embodiment 
of the present subject matter. At 455, an atrial series of characteristic points are 

30 generated from an atrial channel and a ventricular series of characteristic points 
are generated from a ventricular channel. The total number of atrial 
characteristic points is denoted as Ac nt and the total number of ventricular 
characteristic points is denoted as V cn t. At 460, a time window is established in 
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advance of each characteristic point in the ventricular channel. At 465, a i 
backward count is generated by counting, in a single pass through the ventricular 
series of characteristic points, all ventricular characteristic points that have at 
least one atrial characteristic point in the window. At 470, a forward count is 
5 generated by counting, in a single pass through the ventricular series of 
characteristic points, all atrial characteristic points that fall into at least one 
window associated with a ventricular characteristic point. In one embodiment, 
these counts may be performed in the same single pass through the ventricular 
characteristic points. At 475, X and Y values are generated based on the forward 

r 

10 count, the backward count, Acnt and Vent In one embodiment, the X- value 
indicates the backward count expressed as a percentage of all ventricular 
characteristic points. In one embodiment, the Y-value indicates the ratio of 
forward count to backward count. In one embodiment, the Y-value indicates 
(A:V) * Vcnt/Acnt 2 . At 480, the X-Y point is classified to discern a rhythm type. ■ 

1 5 Classification of the data, in one embodiment, includes selecting a separation 
line or separation contour in the X-Y plane. Data falling on one side of the 
separation line or contour denotes a particular rhythm type while data on another 
side denotes a different rhythm type. In one embodiment, a separation line or 
contour is selected according to methods known as Support Vector Method 

20 (SVM) or support vector classification (SVC) (see V. Vapnik. The Nature of 
Statistical Learning Theory. Springer- Verlag, New York, 1995, and others). 
Factors to consider in placement of a separation line include determining the 
placement of the line and a slope. Considerations include the proximity of the 
data points and the a width of a separation region between different portions of 

25 the X-Y plane. In one embodiment, the X-Y plane is regionalized in a manner 
that maximizes the separation of data points. 

Other embodiments are also contemplated. For example, in one 
embodiment, a window is established after each atrial characteristic point and 
counts are established based on the relative timing of atrial characteristic points 

30 and ventricular characteristic points. 

Other metrics are also contemplated for discriminating among rhythm 
types. For example, in one embodiment, discrimination is based on the 
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percentage of atrial characteristic points that are not associated with a ventricular 
characteristic points or a value derived from ventricular characteristic points. 

Exemplary Device 

5 In one embodiment, the present subject matter includes circuitry, 

hardware and software for implantation in a body. In one embodiment, the 
present subject matter is adapted for use external to a body. 

Fig. 14 includes a diagram of one embodiment of portions of cardiac 
rhythm management system 500 according to the present subject matter. As 

10 illustrated, system 500 includes unit 503, lead 510 and lead 540, each of which 
are coupled to heart 505. In various embodiments, unit 503 includes an 
implantable unit and operates as a pacemaker, a CRT device, a 
cardioverter/defibrillator, a pacer/defibrillator and a drug delivery device. In the 
figure, unit 503 senses cardiac activities from heart 505 and delivers therapy to 

1 5 heart 505 according to programming and circuitry of unit 503. Leads 510 and 
540 provide electrical connections between unit 503 and heart 505 and each 
includes an electrode adapted to be disposed in or about heart 505. Lead 510 is 
coupled to input circuit 515 and provides an electrical signal sensed from cardiac 
activity of heart 505. Lead 540 provides therapy to heart 505 as a function of a 

20 signal received from therapy circuit 535. 

In one embodiment, unit 503 includes input circuit 515, controller 520, 
telemetry circuit 525, memory 530 and therapy circuit 535. Input circuit 515 is 
coupled to lead 510 and receives an analog signal based on cardiac activity of 
heart 505. In one embodiment, input circuit 515 includes an analog-to-digital 

25 converter. 

In one embodiment, controller 520 includes a processor and 
programming to implement a method as described herein. In one embodiment, 
controller 520 includes a circuit to control therapy circuit 535 and telemetry 
circuit 525. 

30 Memory 530 is coupled to controller 520 and provides data storage. In 

one embodiment, the input signal is sampled and characteristic points are 
generated substantially in real-time. In one embodiment, characteristic points 
are generated based on stored sample data derived from the input signal. In one 
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embodiment, programming executing on a processor of controller 520 is stored 
in memory 530. 

In one embodiment, controller 520 determines what therapy or electrical 
stimulation is to be delivered by therapy circuit 535 to heart 505 using lead 540. 
5 In one embodiment, controller 520 determines what data is to be 

communicated using telemetry circuit 525. Telemetry circuit 525, in various 
embodiments is adapted to provide wired or wireless telemetry. 

In one embodiment, controller 520 executes programming to determine 
curvature of an input signal and generates characteristic points. 

10 In one embodiment, controller 520 includes a comparator which 

generates an output based on a comparison of an input signal and a reference 
signal. Exemplary comparisons include determining if the lobe is making a 
positive or negative excursion from a baseline. In one embodiment, an 
exemplary comparison includes determining if the curvature series lobe is 

15 greater or less than a dead-zone threshold. In one embodiment, an exemplary 
comparison includes comparing the curvature signal with a hysteresis value of a 
threshold. 

In the figure, two leads are shown coupled to heart 505. In one 
embodiment, a single lead provides both sensing and therapy for heart 505. 
20 System 500, in one embodiment, is adapted to determine a heart rate 

according to the methods described herein. The heart rate is then stored in 

memory 530. In one embodiment, the heart rate is used to determine appropriate 

» 

therapy or to confirm application of a particular therapy. In one embodiment, a 
first heart rate is determined by a method described herein and a second heart 

25 rate is determined by an auxiliary method and controller 520 reconciles the first 
heart rate and second heart rate to determine a third heart rate. In one 
embodiment, the heart rate is communicated to a remote programmer by 
telemetry circuit 525. 

System 500, in one embodiment, is adapted to classify a rhythm 

30 according to the methods described herein. The rhythm classification is then 
stored in memory 530. In one embodiment, the classification is used to 
determine appropriate therapy or to confirm application of a particular therapy. 
In one embodiment, a first classification is determined by a method described 
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herein and a second classification for the same epoch is determined by an 
auxiliary method and controller 520 reconciles the first classification and second 
classification to determine a third classification. In one embodiment, the 
classification is communicated to a remote programmer by telemetry circuit 525. 

5 

Alternative Embodiments 

Characteristic points can be described in terms of a size corresponding to 
the area of a lobe in the curvature series. In one embodiment, before calculating 
a heart rate or classifying a rhythm type according to a method described herein, 

1 0 the series of characteristic points are filtered based on their size. In one 

embodiment, this entails selecting only those characteristic points having a size 
greater than a predetermined value for use in determining a heart rate or rhythm 
classification. In one embodiment, other filtering criteria is applied to the series 
of characteristic points before determining a rate or rhythm. 

15 In one embodiment, a rhythm is classified based on a straight line in the 

X-Y plane. In one embodiment, an irregular shaped region is described in the x- 
y plane and the cluster of data points in the plane relative to the region provides 
the basis for discrimination. In one embodiment, probability functions are 

« 

calculated for the data points and serve as discrimination criteria. 
20 In one embodiment, other discrimination techniques are used in 

conjunction with those described herein and a processor resolves any differences 

and makes a classification selection. 

In one embodiment, a method described herein is executed by hardwired 

circuitry. The circuitry may be analog or digital and include circuit elements 
25 such as operational amplifiers or logic gates. 

Conclusion 

The above description is intended to be illustrative, and not restrictive. 
Many other embodiments will be apparent to those of skill in the art upon 
3 0 reviewing the above description. 



33 



WO 2005/002669 PCT/US2004/020723 

What is claimed is: 

1 . An implantable cardiac rhythm management device comprising: 

an input circuit for receiving a sampled signal corresponding to cardiac 
electrical activity; 

a controller coupled to the input circuit and adapted to determine a 

i 
j 

fundamental frequency of the sampled signal by autocorrelating a function of a 
series of characteristic points, the series of characteristic points determined based on 
the sampled signal and each characteristic point having a time of a lobe in a 
curvature series; and 

a memory coupled to the controller and adapted to store the fundamental ' 
frequency. 

> 

2. The device of claim 1 further wherein the controller is adapted to determine a 
size of a characteristic point and wherein the controller is adapted to determine the 
fundamental frequency as a function of the size. 

3. The device of claim 1 further comprising a rate estimator coupled to the 
controller and adapted to generate a beat frequency and further wherein the 
controller is adapted to generate a heart rate based on the beat frequency and the 
fundamental frequency. 

4. The device of claim 1 further including a telemetry circuit coupled to the 
controller and adapted to communicate to a programmer. 

5. The device of claim 1 further including a therapy circuit coupled to the 
controller and adapted to deliver electrical stimulation as a fiinction of a signal 
received from the controller. 
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6. A method comprising: 

from a curvature series generated as a function of a sampled input signal, 
establishing a series of characteristic points with each characteristic point 
corresponding to a time of occurrence of a lobe in the curvature series; 

using a processor to determine a frequency for the input signal by 
autocorrelating a function of the series of characteristic points; and 

storing the frequency in a memory. 

7. The method of claim 6 wherein the time of occurrence of the lobe includes 
the time of occurrence of a centroid of the lobe. 



8. The method of claim 6 wherein, for each characteristic point, determining a 
size as a function of an area of the lobe, and further wherein using the processor to 
determine the frequency for the input signal includes determining the frequency as a 
function of the size of each characteristic point. 

■ 

9. The method of claim 6 wherein autocorrelating the function of the series of 
characteristic points includes autocorrelating a series of time differences for adjacent 
characteristic points as a function of time of occurrence for each characteristic point. 

■ 

1 0. The method of claim 8 wherein autocorrelating the series of time differences 
for adjacent characteristic points as the function of time of occurrence for each 
characteristic point includes evaluating a product of a time difference between 
adjacent characteristic points and a time difference between time shifted adjacent 
characteristic points. 

1 1 . The method of claim 10 wherein the adjacent characteristic points are in time 
overlap relation with the time shifted adjacent characteristic points. 

12. The method of claim 6 wherein autocorrelating the function of the series of 
characteristic points includes autocorrelating a product of at least two factors 
selected from the series of characteristic points. 
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13. A method comprising: 

receiving a sampled input signal as a function of sensed cardiac electrical 
activity; 

using a processor to generate a curvature series as a function of each received 
sample; 

generating a series of characteristic points as a function of the curvature 
series, each characteristic point associated with a lobe in the curvature series, 
wherein each characteristic point has a time as a function of a time of occurrence of 
the lobe and. a size as a function of an area of the lobe; 

autocorrelating a function based on the series of characteristic points to 
determine a fundamental frequency; and 

storing the fundamental frequency in a memory coupled to the processor. 

1 4. The method of claim 1 3 wherein autocorrelating the function includes 
autocorrelating in a time domain. 

1 5 . The method of claim 1 3 wherein autocorrelating the function includes 
autocorrelating in a characteristic point domain. 

1 6. The method of claim 13 wherein receiving a sampled signal includes 
receiving a ventricular rate electrogram and further wherein autocorrelating the 
function based on the series of characteristic points includes time domain 
autocorrelating a characteristic point time difference. 

17. The method of claim 13 wherein receiving a sampled signal includes 
receiving a defibrillation channel electrogram and further wherein autocorrelating 
the function based on the series of characteristic points includes time domain 
autocorrelating a characteristic point time difference function. 

1 8 . The method of claim 1 3 further including delivering therapy as a function of 
the fundamental frequency. 
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19. A system comprising: 

means for receiving a sampled input signal from cardiac electrical activity; 
means for generating a curvature series based on the input signal; 
means for identifying a series of lobes in the curvature series; 
means for determining a time of occurrence of a centroid for each lobe of the 
curvature series; 

means for establishing a series of characteristic points, with each 
characteristic point corresponding to a lobe; and 

means for determining a fundamental frequency for the input signal by 
autocorrelating a function of the series of characteristic points. 

20. The system of claim 19 wherein the means for receiving includes a cardiac 
electrode. 

2 1 . The system of claim 1 9 wherein the means for generating the curvature 
series includes a processor. 

22. The system of claim 19 wherein the means for determining the time of 
occurrence includes means for calculating an area for each lobe. 

23 . The system of claim 1 9 further including: ' 

secondary means for determining a second frequency as a function of the 
input signal; and 

reconciliation means for generating a beat frequency as a function of the 
fundamental frequency and the second frequency. 

24. The system of claim 1 9 wherein the means for determining the fundamental 
frequency for the input signal by autocorrelating the function of the series of 
characteristic points includes means for determining the fundamental frequency for 
the input signal by autocorrelating a series of time differences for adjacent 
characteristic points as a function of time of occurrence for each characteristic point. 
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25. An article comprising a machine-accessible medium having associated data 
wherein the data, when accessed, results in a machine performing: 
receiving a sampled signal; 

generating a curvature series based on the sampled signal; 

generating a series of characteristic points in the sampled signal, each 
characteristic point corresponding to a lobe in the curvature series and having a time 
corresponding to a time of occurrence of the lobe; 

determining a frequency by autocorrelating a function of the series of 
characteristic points; and 

storing the frequency in a memory. 

.' 26. The article of claim 25 wherein the data, when accessed, further results in the 
machine generating, for each characteristic point, a size determined as a function of 
the area of the lobe in the curvature series. 

27. The article of claim 26 wherein autocorrelating the function of the series of 
characteristic points includes autocorrelating the function of the series of 
characteristic points having the size greater than a predetermined value. 

28. An implantable cardiac rhythm management device comprising: 

a ventricular channel input circuit for receiving a sampled ventricular input 
signal for an epoch; 

an atrial channel input circuit for receiving a sampled atrial input signal for 
the epoch; 

a controller coupled to the ventricular channel input circuit and coupled to 
the atrial channel input circuit, the controller adapted to determine relative timing 
between ventricular characteristic points and atrial characteristic points, each 
characteristic point determined based on a curvature series generated as a function of 
the input signals, each characteristic point having a time of a lobe in the curvature 
series; and 
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a discriminator circuit coupled to the controller and adapted to classify the 
epoch based on a comparison of a first function of the atrial characteristic points and 
a second function of the ventricular characteristic points. 

29. The device of claim 28 wherein the discriminator circuit is adapted to 
generate a plot of the first function of the atrial characteristic points and the second 
function of the ventricular characteristic points. 

30. The device of claim 29 wherein the discriminator circuit is adapted to 
establish a separating contour in the plot. 

3 1 . The device of claim 29 wherein the discriminator circuit is adapted to 
classify the epoch by establishing a separating line in the plot. 

32. The device of claim 3 1 wherein the discriminator circuit is adapted to 
establish a separating region proximate the separating line. 

33 . The device of claim 28 wherein the controller is adapted to determine 
relative timing by generating a window prior to each ventricular characteristic point. 

34. The device of claim 30 wherein a dimension of the window is remotely 
selectable. 

35. The device of claim 30 wherein a dimension of the window is varied 
according to a measured parameter. 

36. The device of claim 28 wherein the discriminator circuit is adapted to 
determine a backward count based on a number of ventricular characteristic points 
having at least one atrial characteristic point in the window and adapted to determine 
a forward count based on a number of atrial characteristic points located within at 
least one window of a ventricular characteristic point. 
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37. The device of claim 36 wherein each characteristic point includes a size 
determined as a function of an area of the lobe, and wherein the discriminator circuit 
is adapted to classify the epoch as a function of the size of each characteristic point. 

38. The device of claim 28 further including a memory coupled to the controller 
and adapted to store the classification of the epoch. { 

39. The device of claim 28 further including a telemetry circuit coupled to the 
controller and adapted to communicate with a programmer. 

40. The device of claim 28 further including a therapy circuit coupled to the 
controller and adapted to deliver electrical stimulation as a function of a signal 
received from the controller. 

41 . A method comprising: 

from a first curvature series generated as a function of a sampled first input 
signal of a cardiac signal, establishing a first series of first characteristic points and 
from a second curvature series generated as a function of a sampled second input 
signal of the cardiac signal, establishing a second series of second characteristic 
points, with each characteristic point corresponding to a time of occurrence of a lobe 
in the curvature series; 

using a processor to determine timewise proximity of the second 
characteristic points relative to the first characteristic points; 

as a function of the timewise proximity, classifying the cardiac signal; and 
storing the classification. 

42. The method of claim 41 wherein the time of occurrence of the lobe includes 
the time of occurrence of a centroid of the lobe. 

43. The method of claim 41 wherein the first input signal is received from a first 
heart chamber and the second input signal is received from a second heart chamber. 
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44. The method of claim 41 wherein the first input signal includes a ventricular 
input signal and the second input signal includes an atrial input signal. 

45. The method of claim 41 wherein classifying includes establishing a plurality 
of windows, with each window located before each first characteristic point. 

46. The method of claim 45 wherein classifying the cardiac signal includes: 

determining a backward count as a function of the number of first 

characteristic points having at least one second characteristic point in a window of 
the plurality of windows; 

determining a forward count as a function of the number of second 
characteristic points located within at least one window of the plurality of windows; 
and 

plotting the backward count as a function of the forward count. 

47. The method of claim 46 wherein plotting the backward count as a function of 
the forward count includes plotting the backward count divided by a total number of 
first characteristic points. 

48. The method of claim 46 wherein plotting the backward count as a function of 
the forward count includes plotting the forward count divided by a ratio of a number 
of first characteristic points in the first series to a number of second characteristic 
points in the second series. 

49. The method of claim 46 wherein plotting the backward count as a function of 
the forward count includes plotting the forward count expressed as a percentage. 

50. A system comprising: 

means for receiving a first sampled input signal and a second sampled input 
signal, each a function of cardiac electrical activity for an epoch; 

means for generating a first curvature series based on the first sampled input 
signal and a second curvature series based on the second sampled input signal; 
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means for identifying a series of lobes in each curvature series; 
means for determining a time of occurrence of a lobe in each curvature 

series; 

means for establishing a first series of characteristic points based on the first 
curvature series and a second series of characteristic points based on the second 
curvature series, each characteristic point corresponding to a lobe; 

means for correlating a time of the occurrence of first series characteristic 
points with a time of the occurrence of second series characteristic points; and 

means for classifying the cardiac electrical activity based on the correlation. 



5 1 . The method of claim 50 wherein the means for determining the time of 
occurrence of the lobe includes means for determining the time of occurrence of a 
centroid of the lobe, 

52. The system of claim 50 wherein the means for determining the time of 
occurrence includes means for determining a size of a centroid of the lobe wherein 
the size is a function of an area of the lobe. 



53. The system of claim 50 wherein the means for correlating includes means for 
establishing a plurality of time windows, each window disposed ahead of a first 
series characteristic point. 

54. The system of claim 53 wherein the means for correlating includes means for 
determining a backward count based on a number of first series characteristic points 
having at least one second series characteristic point within a window of the 
plurality of windows and further includes means for determining a forward count 
based on a number of second series characteristic points disposed within a window 
of at least one first series characteristic point. 

55. The system of claim 54 wherein the means for classifying includes means for 
plotting a function of the backward count and the forward count. 
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56. The system of claim 55 wherein the means for classifying includes means for 
comparing a distribution relative to a separation contour in a plot of the function of 
backward count and the forward count. 



57. The system of claim 50 wherein the means for receiving includes a cardiac 
electrode. 



58. An article comprising a machine-accessible medium having associated data 
wherein the data, when accessed, results in a machine performing: 

receiving a first sampled signal and a second sampled signal, each based on 
cardiac electrical activity for an epoch; 

generating a first curvature series and a second curvature series based on the 
sampled signal; 

generating a first series of characteristic points in the first sampled signal and 
a second series of characteristic points in the second sampled signal, each 
characteristic point corresponding to a lobe in a curvature series and having a time 
corresponding to a time of occurrence of the lobe; 

generating a classification of the epoch based on a plot of timewise 
occurrence of first series characteristic points relative to timewise occurrence of 
second series characteristic points and a separation contour; and 

storing the classification in a memory. 

59. The article of claim 5 8 wherein the data, when accessed, further results in the 
machine generating, for each characteristic point, a size determined as a function of 
the area of the lobe in the curvature series. 

60. The article of claim 58 wherein the data, when accessed, further results in the 
machine generating a plurality of windows, each window disposed ahead of a 
characteristic point of the first series of characteristic points. 

6 1 . The article of claim 60 wherein generating the classification includes plotting 
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■ 

a backward count based on a number of first series characteristic points having at 
least one second series characteristic point within a window of the plurality of 
windows and further includes means for determining a forward count based on a 

* 

number of second series characteristic points disposed within a window of at least 
one first series characteristic point. 
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